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ABSTRACT 

Near the edges of the instability strip the rate of stellar evolution is larger than the growth-rate of the 
pulsation amplitude, and the same holds whenever the star is engaged in pulsational mode switching. 
Stellar evolution therefore controls both the onset of pulsation at the edges of the instability strip and 
of mode switching inside it. Two types of switchings (bifurcations) occur. In a soft bifurcation the 
switching time scale is the inverse harmonic mean of the pulsational modal growth-rate and of the stellar 
evolution rate. In a hard bifurcation the switching times can be substantially longer than the thermal 
time scale which is typically of the order of a hundred periods for Cepheids and RR Lyrae. We discuss 
some of the observational consequences, in particular the paucity of low amplitude pulsators at the edges 
of the instability strip. 
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1. INTRODUGTION 

The Cepheid and RR Lyrae instability strips (ISs) are 
observed to be rich in behavior, as seen for example in 
the beautiful work of Udalski et al. (1999), with stars 
undergoing single mode (SM) pulsations in the funda- 
mental (F) mode, in the first (01) or second overtone 
(02), or beat pulsations in two modes (also called double- 
mode pulsations, DM). Theoretical work has shown that 
for the Galactic Cepheids we have a good understanding 
of the F and 01 Cepheids (Moskalik, Buchler, Marom 
1992, Bono, Marconi & Stellingwerf 2000, Fcuchtinger, 
Buchler & KoUath, 2000) and of RR Lyrae (Feuchtinger 
1999). Furthermore the problem of DM Cepheids and RR 
Lyrae has finally been solved (KoUath, Beaulieu, Buchler 
& Yecko 1998, Feuchtinger 1998, KoUath, Buchler, Szabo 
& Csubry, 2002, hereafter KBSC, KoUath & Buchler, 2001, 
hereafter KB). However some small problems remain, in 
RR Lyrae (KoUath, Buchler & Feuchtinger 2000) and 
Cepheids, in particular those with low metallicity {e.g., 
Feuchtinger et al. 2002). 

As the stars navigate through the IS their structure un- 
dergoes changes. Accordingly, the allowed stable pulsa- 
tional states change on an evolutionary time scale, and 
the stars need to switch their pulsational behavior. Dur- 
ing some stages more than one type of pulsational behav- 
ior is possible. Some of these pulsational states may be 
reachable, while other may not be (because a finite am- 
plitude kick would be required). Stars can therefore be 
in a different pulsational state at the same place on the 
blueward and the redward tracks (hysteresis), as for ex- 
ample in the case of Cepheids and RR Lyrae. Some 
classical Cepheids and RR Lyrae have been observed to 
change their pulsational status on human time scales. For 
example, the pulsation amplitude of Polaris has decreased 
dramatically over several decades (e.g., Kamper & Fernie 
1998). Changes have also occurred in RR Lyrae variables 
{e.g., Clement & Goranskij 1999). The theoretical aspects 



of these mode switchings have been addressed in the past 
{e.g., Stellingwerf 1975, Bono, Castellani & Stellingwerf 
1995), but many of the questions that observations pose 
have remained unanswered. For example, observations in- 
dicate a paucity of low amplitude s-Cepheids fundamental 
Cepheids {e.g., Fernie et al. 1995 for the Galaxy, Beaulieu, 
private comm. for LMC and SMC), and no stars with low 
component amplitudes in the DM region (Beaulieu, priv. 
comm., see figure 4 in Buchler 1998). A similar situation 
occurs for RR Lyrae stars {e.g.. Walker & Nemec 1996, 
Walker 1998, Udalski et al. 1997). These deficiencies of 
small amplitude pulsators seem too large to be solely due 
to observational bias. 

Here we address the problem of mode switching in 
a novel fashion, by taking temporal evolution explicitly 
into account. Hydrodynamical modelling, to be feasible, 
has typically decoupled pulsation from evolution (usually 
called a quasi-static approximation). The decoupling al- 
lowed the calculation of all the possible full amplitude and 
steady pulsational states of a given model at any point 
along the evolutionary track, i.e., of a sequence of models. 
However, this quasi-static approximation can break down 
near the bifurcation points where the star has to switch 
from one behavior to another {e.g., Lebovitz 1990). (A bi- 
furcation point is where there occurs a qualitative change 
in the modal behavior). The reason is that, right at a bi- 
furcation point, the growth-rate of a mode or limit cycle 
can vanish and remain small in the vicinity. This causes 
a temporal delay for the pulsation to achieve the full am- 
plitude that the hydrodynamical modelling would suggest 
- the pulsational behavior is then controlled by evolution. 
In the following we make use of the amplitude equation 
formalism (Buchler & Goupil 1984, Buchler 1993) which 
is perfectly tailored to the study of what is happening in 
the vicinity of the modal bifurcation points, and we show 
that there are observational consequences. 
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2. ANALYTICAL CONSIDERATIONS 

We will be led to distinguish two different types of mode 
switchings that we label soft and hard. A soft bifurcation 
occurs when the nascent mode starts off with infinitesimal 
amplitude and the growth-rate is also infinitesimal. At a 
hard bifurcation the stable fixed point disappears (the cur- 
rent pulsational behavior is no longer allowed) and the star 
is left looking for another stable fixed point {i.e., another 
pulsational behavior). The rate of change of the amplitude 
is generally not infinitesimally small in this case. In §3 we 
will see examples of the two types of switchings. 

2.1. Soft Bifurcations 

2.1.1. Hopf Bifurcation 

We first consider the simplest soft modal switching 
which occurs when the star enters or exits the IS. With- 
out lack of generality let us consider the neighborhood of 
the blue edge where the linear growth-rate k of a mode, 
the fundamental or the first or second overtone depending 
on the luminosity, changes from negative (vibrationally 
stable) to positive (unstable). Generically, in the vicin- 
ity of the blue edge K{t) varies proportionately to time, 
K = x^i where the origin of time is taken to be the pas- 
sage through the blue edge. The quantity x = dn/dt is 
therefore the stellar evolutionary rate of change of the lin- 
ear growth-rate. This type of bifurcation is well known 
to dynamicists and goes under the name of supercritical 
Hopf bifurcation {e.g., Berge et al. 1984). The behavior 
of the amplitude is governed by an amplitude equation of 
the form {e.g., Buchler & Goupil 1984, CouUet & Spiegel 
1984) 



j^A^K{t)A-q{t)A^ 



ith 



n{t)=xt (1) 



where K{t) is the instantaneous linear growth-rate of the 
mode and q{t) a nonlinear saturation term which varies 
little along an evolutionary track. 

If the star were perturbed away from its quasi-static 
stable limit cycle amplitude 



A^,{t) = ^Jn{t)lq{t) 



(2) 



(as obtained from setting dA/dt = in Eq. |l|), it would 
return to it on a tth^ time scale, i.e., over a few 
hundred pulsations, because the relative growth-rates are 
typically of the order of a percent inside the IS. Away from 
the bifurcation, the thermal time scale tth is much shorter 
than the stellar evolution time scale tev ~ {dlnTcs/dt)~^ . 

Right at the bifurcation, however, k vanishes by defini- 
tion, and it remains small in some neighborhood, so that 
tth ^ tev for a while, and the amplitude does not grow. In 
fact the star must wait for evolution to increase the value 
of K. Only then can the amplitude catch up with its quasi- 
static value Aqs{t). It should be clear from this discussion 
that the evolutionary rate of change of the growth-rate x 
will play a crucial role. 

Let us call t the time when the solution catches up, and 
A the amplitude at that instant. We are interested in de- 
termining the dependence of t and A = Aqs{t) on x, q and 
Ao=A{to). The latter is the pulsation amplitude at to that 
arises because of noise in the star. 



The coefficient q is time-dependent because of changes 
in the structure of the star as it evolves, but its time- 
dependence is small. In a later example of an actual RR 
Lyrae model we include its time-dependence. Here, we ap- 
proximate q 'AS a constant. With the substitution E ^ A^ 
we can then readily solve Eq. (|^) in the form 



/Eo + 2q / e-^(*'-*")di' 



e-^' /Eo + q/{xt) 



(3) 



where the last expression is the first term of an asymptotic 
expansion for the integral for large x^^ ■ 

The point at which the pulsation amplitude catches up 
with the quasi-static value can be obtained to a good ap- 
proximation by requiring that the first term on the rhs of 
Eq. (||), which completely dominates at first, becomes a 
fraction ^ of the second term (^ = 0.1 turns out to be a 
good value) . With the help of Eq. (^ this yields a simple 
transcendental equation 



2u — In u — In / = , where u = xt , 



and / = 



X 



{iqAlf 



(4) 
(5) 



A good approximation to the solution, when f is large, is 
obtained with the fitting expression 



u = xP ^ 1.2 Log/ + 0.7. 



(6) 



Finally we have t — {u/x)^- The quantity q can be 
re-expressed in terms of the maximum pulsation ampli- 
tude in the IS, since qA^^,^. = Kmax- Thus A/Amax = 

{Xi/l^max)^. 
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Fig. 1; Temporal behavior of the pulsation amplitude after the 
star enters the IS (solid line), for the case of x = 10~^^s~^. 
Curves are for three values of initial perturbation (noise level) 
Logy4o — —14, —16, —18. The dotted line represents the quasi- 
static amplitude Aqs {t) that would be achieved if the pulsation 
could react infinitely fast to evolutionary changes in k. The 
amplitude scale is set by the choice q — 10~^. The thick dots 
correspond to the A obtained with the approximation of Eq. 6. 
The dashed line represents exp(^x^^) /C for the middle value 
of Ao. 

The typical behavior of A{t) as computed from Eq. (|^) 
is plotted as solid lines in Fig. 1 for Logx = — 16 and for 
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three values of Ao (from left to right: Logi^o = -14, - 
16, -18). The quantity q which scales the amplitude, has 
been arbitrarily set to q = 10~^. Shown as a dotted line 
is the quasi-static limit cycle amplitude Aqs{t) — \P(tlq^ 
namely the value that the amplitude would achieve if the 
system could react instantly. 

We have no knowledge of the intensity of the fluctuation 
Ao in a star which represents the projection of the turbu- 
lent noise on the eigenvector of the mode. As suggested by 
the approximation formula, the larger is, the faster the 
amplitude grows (Fig. 1). But one notes that the effect 
is only logarithmic and therefore a precise knowledge of 
Aq is of little importance. It is the rate of change of the 
growth-rate, x, that sets the switching time scale, i.e., the 
scale of the abscissa in Fig. 1. 

The three fat dots denote the t and A for the three val- 
ues of Aa as given by the fitting formula. It can be seen 
that they provide an excellent agreement with the numer- 
ical integrations of Eq. (|I]). The dashed line represents 
Aq exp(ixf^)/^ for the middle value of Aq. 

2.1.2. SM to DM Bifurcation 

When two (nonresonant) excited modes are involved 
Eq. needs to be replaced by the two coupled ampli- 
tude equations for modes 1 and 2 {e.g., Buchler & Goupil 
1984) 

Ai = (^Ki{t) + qiiit)Al + qi2{t)Al + si{t)AlAl^ Ai (7a) 

A2 = (k2(0 + q2l{t)Al + q22(t)Al + S2(t)AlAtj A2 (76) 

We have kept some of the quintic terms in these equa- 
tions which are necessary in order to describe the observed 
modal behavior (KB, KBSC). (Earlier studies by Buchler 
& Kovacs (1987) had hoped to get away with retaining 
only the cubic ones). 

When the transition from SM to DM inside the IS is 
such that the nascent mode starts with infinitesimal am- 
plitude, as we assume in this section on soft bifurcations, 
then the behavior of this amplitude is very similar to the 
Hopf bifurcation that we have just described. The reason 
is easy to see: in the vicinity of the bifurcation the ampli- 
tude of the finite amplitude mode (say Ai) changes very 
little, and we can expand around this point which leads to 
an equation for the amplitude say A2 of the nascent mode 
that is identical to Eq. ^ 

To summarize then, stellar evolution causes a delay of 
order during which the pulsation amplitude stays neg- 
ligibly small, but after which it very rapidly achieves its 
quasi-static value. The time spent at the small interme- 
diate amplitudes is therefore very small and the stars are 
very unlikely to be observed with such small amplitudes. 
The amplitude at catch-up A^^ x^ thus also depends on 
and when the stellar evolution is slow it will be large. Fi- 
nally, we note that the time scale and the catch-up ampli- 
tude are both insensitive to the poorly known noise level, 
typified by A^. 

2.2. Hard Bifurcations 

A hard bifurcation occurs when, as a result of stellar 
evolution, a fixed point merges with another unstable fixed 



point, in a mutual annihilation. One can readily visual- 
ize this in an amplitude-amplitude plane (see e.g., KB, 
KBSC), where the fixed points, stable and unstable, lie at 
the intersection of two loci (curves). As Tcff changes with 
evolution, the two intersection points coalesce and then 
disappear as the curves move apart. At that point the star 
has no choice but to look for another stable pulsational 
state to mew into. The new type of pulsation depends on 
the overall topography at the time, and is determined by 
which fixed point attracts the stellar model. One might 
think that the time to attain the new pulsational state is 
the thermal time scale tth — n^^, but it can be and it 
is generally much longer as we shall show in a realistic 
example. We note that for a hard bifurcation the stellar 
evolution time scale plays no important role because the 
amplitude growth-rates do not vanish and thus the time 
scale for amplitude change is always finite (KBSC). 

3. CEPHEIDS AND RR LYRAE 
3.1. The Blue Edge 

Let us consider first a redward evolving 6 Mq Cepheid 
that enters the first overtone IS. We use Eq. (^) to estimate 
the minimum observable amplitude A. Stellar evolution 
calculations indicate that a 6Mq Cepheid variable evolves 
at a rate of dT^s / dt ^ 2xlQ~^'^ K when entering the 
overtone blue edge {e.g., Alibert et ai 1999). Linear mod- 
els give a rate of change dn/dTcs^ 2xlO~^'^{Ks)^^ , so 
that we obtain % ~ 4x 10~^°s^^. The maximum growth- 
rate Kmax in the IS is 4xlO^*s^^. If we assume intrin- 
sic fluctuations of the order of Ao/Amax^iO~^'^, where 
Amax is the maximum amplitude in the instability strip, 
we find that the u in Eq. 6 takes on values around 35 - 
40. Thus we estimate that A/A„iax ^0.1 at the blue edge, 
and that the pulsation amplitude stays negligibly small 
for t ^ 1000 yr, and that it then increases very rapidly to a 
value 10 % of the maximum amplitude in the IS. Because 
of this rapid increase of the amplitude very few overtone 
Cepheids with amplitudes less than A ^ Amax should be 
observed near the blue edge. 

Turning now to a 8Mq Cepheid which evolves at 
dT^sldt^2xlQ-'^Ks-^ (Alibert et al. 1999) when en- 
tering the fundamental blue edge, with a dn/dTcS''^ 
4xl0""(/i:s)-i we obtain x~8xl0-2°s~2. The maxi- 
mum Kmax in the IS is ~ 7x 10~^s~"'^. Thus we estimate 
that again A/ Amax ^O-l near the blue edge, and that the 
pulsation amplitude stays negligibly small for i~800yr. 

Cepheids evolve both blueward and redward. Only that 
half of the Cepheids near the blue edge will therefore be 
affected, namely only those that move redward and enter 
the IS. For the same reasons, there should also be a (50%) 
deficiency of Cepheids at the red edge. Furthermore at 
a given period there will be a medley of Cepheids spread 
over over the whole IS. Consequently there should be an 
overall deficiency of low amplitude Cepheids at all periods. 

For RR Lyrae stars we deduce from Dorman (1992) 
that dTcs/dt^ 1.2xl0~^^s~^ when entering the overtone 
blue edge. From our linear models we obtain dn/dTcs 
l.5xlO-^{Ksy^ which leads to x ~ l-8x 10~2is"2_ 
the IS Kmax is 71 X 10~^s~^. We thus estimate that 
A/Amax "^0-0^ at the blue edge, and that the pulsation 
amplitude stays negligibly small for i ~ 5000 yr. 



3.2. Switching to and From Double-Mode Pulsation 

The coefficients of the amplitude equations (|^) can be 
extracted from hydrodynamical simulations, as we will il- 
lustrate for a sequence of RR Lyrae models. First we com- 
pute the stable pulsational states at selected values in the 
Tcff range of interest. From these we extract for each model 
(characterized by Toff) the coefficients Ko(Toff), goo(7oCf), 
.... The values of d In Teff /dt along the evolutionary track 
through the IS determine the time-dependence of Ko{t), 
qoo{t), .... 
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Fig. 2: Redward evolution of the RR Lyrae model through the 
HR diagram. Top: fundamental (F) SM amplitude or com- 
ponent amplitude in the DM; bottom: Ol SM amplitude or 
component amplitude. At (a): switching from SM Ol to DM 
pulsation and at (/3): switching from DM to the F limit cycle. 

We illustrate the method for an RR Lyrae star as it 
evolves through the switching from SM overtone to DM to 
SM fundamental pulsation. Figure 2 shows the results of 
the integration of Eqs. (0), i.e., the pulsation amplitudes 
of the F and 01 modes. For convenience we have cho- 
sen a track with fixed luminosity h— 50Lq and mass M= 
0.77Mo, Z= 0.0001, rather than an inclined evolutionary 
track (Dorman 1992). The star is assumed to be evolving 
toward colder Tgff. The origin of time is arbitrarily taken 



to be a few thousand years before the point (marked a) 
where the 01 limit cycle becomes unstable to a DM pul- 
sation (point P). 

Fig. 2 is best understood in conjunction with paper 
KBSC. (Their Fig. 3 shows the location of the SM and DM 
fixed points and their stability in an amplitude-amplitude 
plot for selected T^gs, while their Fig. 2 shows the vari- 
ation of the quasi-static amplitudes along the sequence). 
The solid lines here denote the pulsation amplitudes, or 
the component amplitudes in the case of DM pulsations, 
as a function of time, top graph for the F amplitude and 
bottom for the 01 amplitude. The dotted lines denote 
the quasi-static amplitudes, i.e., the instantaneous fixed 
points of Eqs. (0). 



T 1 r 



-8 -6 -4 

Log(dA/dt) [((5R/R.)/d] 




0.02 0.04 
Ao [6R/K] 

Fig. 3: F amplitude - Ol amplitude phase space for two RR 
Lyrae models with a lOK difference in effective temperature, 
where the upper one has the higher temperature. The short 
lines represent the normalized flow field (with the dots repre- 
senting the bottom of the vectors). The large solid/open circles 
in the top graph denote the location of the stable/unstable fixed 
points. The square in the bottom denotes the location of the 
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just vanished pair of DMs. The speed of the flow is represented 
by a gray scale. The amplitudes refer to (SR/Rj^urS- 

3.2.1. Soft Bifurcation 

The first inset in Fig. 2 is a blow-up of the onset of 
DM behavior near point (a). As expected, the behavior of 
the nascent F amplitude is very similar to the solution of 
Eq. (0) that was shown in Fig. 1. The major difference is 
that here x becomes the rate of change of the linear stabil- 
ity coefficient of an F mode perturbation of the 01 limit 
cycle (the Floquet stability exponent) rather than that of 
the linear growth-rate of the mode. 

One sees that, for DM onset as well, stellar evolution 
inhibits the growth of the pulsation amplitude for some 
8000 yr, i.e., much beyond the thermal time scale which is 
of the order of a hundred days. DM RR Lyrae (RRd) with 
small fundamental mode component pulsation amplitudes 
should be rare. 

3.2.2. Hard Bifurcation 

An example of a hard bifurcation occurs in the second 
bifurcation for our RR Lyrae model, marked (3 in Fig. 2. 
The stable DM is annihilated by a second unstable DM 
and they both disappear (c/. Figs. 2 and 3 of KBSC). In 
this case the star gets attracted by the F fixed point, and 
thus the pulsation switches to a stable F limit cycle. 

Figure 3 depicts the nature of the two RR Lyrae stel- 
lar models of Fig. 1, taken on either side of point /?, just 
before (upper panel) and after (lower panel) the DM fixed 
point vanishes. In this amplitude-amplitude diagram the 
large solid circles represent the stable fixed points, of which 
now only the F fixed point is left in the lower figure with 
amplitude Ao= 0.055. The large open circle refers to the 
unstable 01 fixed point at Ai — 0.030. The large square 
denotes the location of the just vanished pair of DMs. 

The 'velocity' flow field in the Aq~Ai plane is given by 
u = (Ao, Ai), and the speed by u = {A^ + Al)i . The 
'arrows' in Fig. 3 represent the normalized 'velocity' flow 
field , i.e., the vectors u/u. The dots denote the bases 
of the vectors. The gray scale indicates the Log of the 
speed (length of the velocity vectors). Note that the typ- 
ical 'thermal' time scale is given by the shading near the 
origin. 

It is very apparent that all the trajectories get attracted 
very rapidly to an arc going from the Ai axis to the the 
^0 axis. This arc is an integral curve of the steady state 
amplitude equations (known as a heteroclinic connection 
because it connects two different fixed points, viz. the O 
and F fixed points.) Along this arc the motion is consider- 
ably slower than thermal. The transition occurs on a time 
scale of 100 yr, very long compared to the thermal time 
scale of ^ 100 d. 

At the moment the DM disappears, the star finds itself 
at the point marked with a large square (bottom graph), 
and it gets slowly attracted to the only stable fixed point, 
namely F. The pulsation slowly changes from DM to fun- 
damental SM. 

In the case of a hard bifurcation, the switching time, 
even though of the 'same order of magnitude' as the ther- 
mal time scale, thus can exceed it by factors of ten. 



3.2.3. Speed Along Heteroclinic Connection 

Finally we examine why the speed is slow along the het- 
eroclinic connections that link the SM and the DM fixed 
points. In Fig. 3 one visualizes the location of these hetero- 
clinic connections from the vector field. (They are shown 
explicitly in Fig. 1 of KBSC.) Consider now the two loci 
that are defined by setting the first and second of eqs. @) 
equal to zero, respectively. On these loci either the hor- 
izontal (Ai) or the vertical {A2) speed vanishes by defi- 
nition. These loci each emanate from one of the two SM 
fixed points, and both of them intersect the heteroclinic 
connection at the two DM points. The heteroclinic con- 
nection lies therefore very close to the two loci, especially 
near and in between the DM fixed points, and the speed 
u — (Ai)^ -I- (^2)^)5 must therefore also remain small. 

One also notes that, when both modes are linearly un- 
stable, which is the situation we are interested in, the het- 
eroclinic connection is always attractive sideways. This 
remains true a fortiori when the loci intersect and give 
rise to a DM fixed point (in which case the speed vanishes 
at the DM). 

Note also that the linear stability root of the F fixed 
point along the Aq axis is given by —2kq and, since kq is 
typically large the speed along Aq is high. Mutatis mu- 
tandis one shows that the speed along the Ai axis is also 
high. Although this only concerns the motion along the 
coordinate axis. Fig. 3 shows that it remains true in be- 
tween. 

4. CONCLUSION 

We have shown that there are two types of bifurcations 
with different characteristics. The first is a soft bifurcation 
where the incipient mode starts with infinitesimal ampli- 
tude whose growth is determined by the evolutionary rate 
of change x = dn/dt of the appropriate linear growth-rate 
K. This time scale x~2 can be very large compared to the 
thermal time scale depending on the stellar evolution 
time scale. The second is a hard bifurcation which occurs 
when the rug is pulled from under the pulsation, as it were, 
i.e., when the pulsation must change, because, due to an 
evolutionary change in the star's structure, it is no longer 
allowed or possible. Here the time scale remains thermal, 
but generally with a large multiplication factor. This delay 
and then abrupt onset of pulsation has as a consequence 
that stellar pulsators with low amplitudes should be rarely 
observed, in agreement with observations. 

Another interest of this paper is that it shows the possi- 
bility of obtaining valuable information about stellar struc- 
ture from the observations of a star that is in the process of 
starting to pulsate or of switching its pulsational state. In 
a subsequent paper we will show how the combined use of 
numerical hydrodynamics, time-series analysis, the ampli- 
tude equation formalism and the knowledge of the stellar 
evolution time scale (dTcff /dt) allow one to compute con- 
sistently the evolution of a pulsating star throughout the 
whole IS. 

This work has been supported by NSF (grant AST 
9819608) and by OTKA (T-038440). We also wish to 
thank an anonymous referee for his suggestions. 
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